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ENERGY  TRANSFER  AND  TRIADIC  INTERACTIONS  IN  COMPRESSIBLE 

TURBULENCE* * * § 


F.  BATAILLE*,  YE  ZHOU* ,  AND  JEAN-PIERRE  BERTOGLIO& 

Abstract.  Using  a  two-point  closure  theory,  the  Eddy-Damped-Quasi-Normal-Markovian  (EDQNM) 
approximation,  we  have  investigated  the  energy  transfer  process  and  triadic  interactions  of  compressible 
turbulence.  In  order  to  analyze  the  compressible  mode  directly,  the  Helmholtz  decomposition  is  used.  The 
following  issues  were  addressed:  (1)  What  is  the  mechanism  of  energy  exchange  between  the  solenoidal  and 
compressible  modes,  and  (2)  Is  there  an  energy  cascade  in  the  compressible  energy  transfer  process?  It  is 
concluded  that  the  compressible  energy  is  transferred  locally  from  the  solenoidal  part  to  the  compressible 
part.  It  is  also  found  that  there  is  an  energy  cascade  of  the  compressible  mode  for  high  turbulent  Mach 
number  ( Mt  >  0.5).  Since  we  assume  that  the  compressibility  is  weak,  the  magnitude  of  the  compressible 
(radiative  or  cascade)  transfer  is  much  smaller  than  that  of  solenoidal  cascade.  These  results  are  further 
confirmed  by  studying  the  triadic  energy  transfer  function,  the  most  fundamental  building  block  of  the  energy 
transfer. 

Key  words.  Compressible  turbulence,  energy  transfer,  triadic  interactions,  closure  theories. 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  It  is  well  known  that  compressible  turbulence  plays  a  prominent  role  in  a  wide  range 
of  important  scientific  and  engineering  applications,  including  high  speed  transport,  supersonic  combustions, 
and  acoustics.  Recently,  a  large  body  of  publications  have  been  devoted  to  study  various  aspects  of  com¬ 
pressible  turbulence  using  both  direct  numerical  simulations  (DNS)  and  large  eddy  simulations  (LES).  Here 
we  simply  mention  a  few  representative  works,  such  as  Feiereisen  et  al.  (1981),  Passot  and  Pouquet  (1985), 
Lee  et  al,  (1991),  Sarkar  et  al.  (1991),  Erlebacher  et  al.  (1992),  Kida  and  Orszag  (1990),  Blaisdell  et  al. 
(1993)  and  Porter  et  al.  (1994).  For  a  comprehensive  review,  the  reader  is  referred  to  Lele  (1995).  These 
numerical  simulations  have  substantially  improved  our  understanding  of  compressible  turbulence.  Never¬ 
theless,  some  basic  physical  processes  of  compressible  turbulence,  such  as  the  energy  transfer  and  triadic 
interactions,  have  not  been  explored  even  at  low  turbulent  Mach  number.  For  example,  do  we  expect  an 
energy  cascade  process  of  compressible  velocity  modes?  How  does  the  energy  exchange  process  between  the 
solenoidal  and  compressible  modes  take  place?  These  type  of  studies  require  a  substantial  spectral  scale 
range  of  interactions.  As  a  result,  it  is  very  hard  to  utilize  the  DNS  databases  since  these  simulations  are 
limited  to  very  low  Reynolds  numbers  and  have  only  very  limited  spectral  ranges.  While  LES  can  provide 
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databases  at  substantially  higher  Reynolds  numbers,  subgrid  models  must  be  introduced. 

Another  way  to  generate  high  Reynolds  number  databases  is  by  using  two-point  closure  models.  Di¬ 
rect  Interaction  Approximation  (DIA)  of  Kraichnan  (1959)  is  a  well-established  approach.  Many  authors 
have  studied  incompressible  turbulence  using  DIA  (see  Leslie  (1973)).  The  method  of  Eddy- Damped- Quasi- 
Normal-Markovian  (EDQNM)  models  (Orszag  (1970))  has  been  shown  as  a  simpler,  but  effective,  alternative 
to  DIA.  Recently,  Bertoglio  et  al.  (1996)  have  presented  DIA  and  EDQNM  equations  for  a  weakly  compress¬ 
ible  turbulence. 

The  aim  of  this  paper  is  to  use  the  EDQNM  closure  theory  to  study  the  energy  transfer  and  triadic 
interactions  of  compressible  turbulence.  The  paper  is  organized  as  follows:  First,  we  review  closure  assump¬ 
tions  of  EDQNM  and  present  the  resulting  transport  equations.  Second,  we  perform  a  detailed  analysis  of 
the  non-linear  transfer  terms.  Finally,  we  investigate  the  most  fundamental  aspect  of  the  energy  transfer 
process,  the  incompressible  and  compressible  triadic  interactions. 


Figure  l:Solenoidal  and  compressible  velocity  correlations 
for  different  turbulent  Mach  numbers 


Fig.  1. 

2.  Compressible  EDQNM  Model.  The  basic  set  of  equations  are  the  Navier-Stokes  and  continuity 
equations.  The  fluid  is  assumed  to  be  homogeneous,  isotropic  and  barotropic.  Reynolds  average  and  Fourier 
transform  are  used  to  obtain  the  fluctuating  turbulent  field  equations  in  spectral  space.  The  equations  are 
partially  linearized  with  respect  to  the  density  fluctuation  which  leads  to  the  condition  Mt  <  1,  where  Mt  is 
the  turbulent  Mach  number  defined  as  Mt  =  yfq^/c o*  Q2  is  twice  the  turbulent  kinetic  energy,  and  cq  is  the 
sound  speed.  To  analyze  the  compressibility  effects,  we  use  the  Helmholtz  decomposition  to  split  the  velocity 
vector  into  a  solenoidal  part  us(K ,  t),  which  corresponds  to  the  velocity  fluctuations  perpendicular  to  the 
wave  vector  K  in  the  Fourier  space,  and  a  compressible  part  uc’(K,  £),  which  corresponds  to  fluctuations  in 
the  direction  of  the  wave  vector. 

The  classical  DIA  approach  (two  point  -  two  times)  is  used  to  derive  compressible  DIA  equations.  These 
equations  are  the  starting  point  used  to  derive  the  compressible  EDQNM  equations  (in  two  points  and  one 
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time).  The  methodology  is  the  same  as  that  of  the  incompressible  case  (Leslie  (1973)  and  Lesieur  (1987)) 
but  more  equations  are  now  needed.  For  more  details,  see  Bertoglio  et  al.  (1996)  or  Bataille  (1994). 

2.1.  Spectral  Equations.  The  EDQNM  governing  equations  of  weakly  compressible  turbulence  are 
the  following: 

-  an  equation  for  the  spectrum  ( Ess ),  the  auto-correlation  of  the  solenoidal  part  of  the  velocity  field: 

(1)  ^ Ess(K,t )  =  -2vK2Ess(K,t)+Tss(K,ty, 

-  an  equation  for  the  spectrum  ( Ecc ),  the  auto-correlation  of  the  “purely  compressible”  part  of  the 
velocity  field: 

n 

(2)  ^  Ecc (K,  t)  =  -2 is'K2Ecc(K,t)  +  Tcc(K,t)  -  2caKEcp(K,t ); 

-  an  equation  for  the  spectrum  of  the  potential  energy  ( Epp )  associated  with  the  pressure: 


(3)  ^  Epp(K,t)=2c0KECP(K,t)- 

-  an  equation  for  the  spectrum  of  the  pressure- velocity  correlation  ( ECP ): 

(4)  ^  ECP(K,t)  =  -iy'K2Ecp(K,t)  +  Tcp(K,t)  +  c0K(Ecc(K,t)  -  Epp(K,t)). 


In  the  case  of  a  Stokesian  fluid: 

(5) 


A  -\~  2(i  4 

<  p>  ~~  3 


p  and  A,  two  dynamic  viscosities,  are  assumed  to  be  uniform. 


2.2.  Energy  Transfer  Terms.  In  equations  (l)-(4),  Tss ,  Tcc,  and  TCP  are  the  transfer  terms.  They 
contain  several  contributions: 


(6) 

(7) 

(8) 


rpSS  _  rpSS  _j_  rpSS  _j_  rpSS  _j_  rpSS  _j_  pS S 


pCC  _  pCC  _|_  pCC  pCC  _|_  pCC  _|_  pCC  _|_ 


tcp  =  Tfp  +  r2CP  +  t3cp  +  rfp  +  r5CP  +  r6cp 

Different  contributions  appearing  in  the  transfer  term  acting  on  the  solenoidal  field  are: 


(9) 


(10) 


T?s 


rpSS 
1  2 


[  & 

Ja  PQ 


1  —  xyz  —  2y2z 


2,2 


f  K3  (1  -y2){x2  +y2) 

Ja  PQ  1  -  x2 


OsksPqS~ssEss(P,  t)Ess  (Q,  t)dPdQ 
0hKPQS~CCEss(P,  t.)Ecc (Q,  t)dPdQ 


(11) 


Tis  =  -  J  ^  (xy  +  z3)esKSp^s~ssEss(K,  i)Ess(Q,  t)dPdQ 
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(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 


T?s  =  J  ^  (2xy)6sKsp-ss-ccEss(K,  t)Ecc(Q,  t)dPdQ 

Tis  =  -  J  ~  (z(l  -z2))0f?f£c-ssEss(K,t)Ess(Q,t)dPdQ 
Different  contributions  to  the  transfer  term  in  the  Ecc  equation  are: 

Tfc(K,t)  =  J  ^  ((x  +  yz)2)8^%ss~ssEss(P,t)Ess(Q,t)dPdQ 

T?c(K,t)  =  jf  ^  "JiZSf  °TpQSS~CCESS(P,t)Ecc(Q,t)dPdQ 

Tfc(K,t)  =  jf  {x2)6ckcp-qcc-ccEcc{P, t)Ecc(Q, t)dPdQ 

Tfc(K,  t)  =  -  J  ^  2z(l  -  z2)ecK%ss-ssEcc{K,  t)Ess(Q,  t)dPdQ 

T°C(K,  t)  =  -J^~  (2z3  -z  +  xy)d%%cc-ssEss(Q,  t)Ecc{K,  t)d,Pd.Q 

TqC(K,  t)  =  f  ^  (2xy)0%%cc-ccEcc(K,t)Ecc(Q,t)dPdQ 
And  lastly,  the  contributions  to  the  transfer  term  in  the  pressure  velocity  correlation  are: 

T?p(K,t)  =/Ag  Opcpqss-ssEss(P  t)Ess(Q,  t)dPdQ 

T?p(K,t)  =  Ja^  6p%ss-ccEss(P, t)Ecc(Q, t)dPdQ 

T?p(K,t)  =  ^  y  8p%cc-ccEcc(P,t)Ecc(Q,t)dPd.Q 

Tfp(K,t )  =  1  -  z2)Bpk%ss-ssEPc{K,  t)Ess(Q,  t)dPdQ 

T°p(K,t)  =  -j T  ^  —  Zl+jH  8p%cc-ssEss(Q,t)Ecc(K,t)dPd,Q 
T§P{K,  i)  =JA^-  {xy)8p%cc-ccEcc{K,t)Ecc(Q,t)dPd.Q 
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Figure  2:Solenoidat  transfer  term 
FlCl.  2. 

The  integration  in  the  P,  Q  plane  extends  over  a  domain  such  that  K,  P  and  Q  form  a  triangle.  The 
expressions  x,y,z  are  standard  coefficients  associated  with  the  geometry  of  the  triad  and  they  the  cosines 
of  the  angles  respectively  opposite  to  K,  P,  Q.  Furthermore,  temporal  functions  (defined  by  0’s)  are  the 
decorrelation  times  deduced  by  integrating  the  DIA  equations  over  time.  Their  expressions  can  be  found  in 
Bertoglio  et  al.  (1996). 

The  transfer  terms  appearing  at  the  transport  equations  for  the  solenoidal  and  compressible  modes 
will  be  analyzed  in  detail  in  the  next  section.  Since  there  is  no  transfer  term  in  the  equation  for  the 
pressure- pressure  correlation,  this  equation  characterizes  the  exchange  between  the  solenoidal  mode  and  the 
compressible  mode  ( Ecc  +  Epp) .  Finally,  the  effect  of  Tcp  is  to  interchange  energy  between  Ecc  and  Epp 
(Bataille  and  Bertoglio  (1993)). 

2.3.  Energy  Spectra.  A  force  is  applied  to  the  large  scales  of  the  solenoidal  velocity.  Our  analysis  is 
carried  out  when  both  solenoidal  and  compressible  modes  have  reached  their  asymptotic  stationary  states. 

In  Figure  1,  we  present  the  spectra  of  both  solenoidal  and  compressible  components  of  the  velocity 
correlation  corresponding  after  the  asymptotic  state  is  reached.  The  Taylor  micro-scale  Reynolds  number, 
Re ,  is  approximately  140.  Here  Rc  =  and  e  is  the  dissipation.  In  the  figures,  the  dimensional  unit  of 

the  correlation  spectra  is  given  as  m?s~ 2  and  the  dimensional  unit  of  wavenumber  K  is  defined  as  m~l.  The 
energy  associated  with  the  purely  compressible  mode  is  found  to  vary  as  the  square  of  the  turbulent  Mach 
number.  At  low  and  moderate  Mt,  the  spectrum  of  the  compressible  component  shows  a  K“n/3  behavior 
in  the  inertial  range  (Bataille  and  Bertoglio  (1993)  and  Bertoglio  et  al.  (1996)).  This  behavior  has  been 
confirmed  with  Large  Eddy  Simulation  by  Bataille  et  al.  (1996). 

In  order  to  examine  the  dependence  of  the  turbulent  Mach  number,  it  was  allowed  to  vary  from  10-4  to 
1.  Note,  however,  strictly  speaking,  our  model  is  valid  only  for  small  Mt.  Nevertheless,  computations  were 
carried  out  up  to  Mt  =  1  in  order  to  study  the  limiting  behavior  of  the  model. 


5 


T 

T  3 

— 0 — 

1  *-m-m  ^ 

— i —  X2 

>1L 

Fig.  3. 


Fig.  4. 


3.  Study  of  the  Transfer  Terms. 

3.1.  Solenoidal  Transfer  Term.  In  Figure  2,  we  observe  that  T'5'5  has  the  usual  shape  observed  in 
incompressible  turbulence  studies.  Specifically,  Tss  is  negative  in  the  large  scales  and  positive  for  small 
scales.  Physically  this  corresponds  to  the  energy  transfer  from  the  large  scales  to  the  smaller  ones. 

There  are  17  contributors  to  the  transfer  function  of  compressible  turbulence.  This  should  be  be  com- 
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figure  5:SolenoidaI  transfer  term  at  different  Mt 


Fig.  5. 


Fig.  6. 


pared  with  its  counterpart  of  incompressible  turbulence: 

(26)  TSSS(K,  t )  =  T?S(K,  t)  +  Tfs(K,  t). 

Indeed,  the  two  key  contributors  to  Tss  are  Tfs  and  T£s ,  which  are  the  same  terms  that  one  finds  in 
incompressible  turbulence.  These  terms  are  usually  called  “input”  and  “output”  terms.  For  our  weakly 
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Fig.  7. 


Fig.  8. 

compressible  turbulence,  they  are  much  more  important  contributors  than  the  others  (see  Figure  3).  Indeed, 
other  contributors  are  negligible  in  comparison  with  Tfs  and  T^s .  We  stress  that  the  summation  of  these  two 
terms  is  the  net  energy  transfer  of  incompressible  turbulence  (solid  line).  To  clearly  illustrate  the  behaviors 
of  the  new  terms,  we  plotted  these  “compressible”  contributors  (Tlf5’,  Tfs  and  T§s)  in  an  enlarged  scale 
(Figure  4)  in  order  to  observe  their  relative  magnitudes.  We  found  that  T«f's  is  much  larger  than  T^s  and 
Tfs,  and  consequently  it  makes  the  largest  contribution  to  the  compressible  effects  in  Tss .  An  important 
feature  of  this  term  is  that  it  is  negative  for  all  spectral  space,  indicating  that  the  energy  transfer  has  been 
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Figure  9:  Compressible  transfer  term 
Fig.  9. 


Fig.  10. 

transferred  from  the  solenoidal  mode  to  the  compressible  mode. 

We  observe  in  Figure  5  that  the  compressibility  has  negligible  effect  on  the  solenoidal  transfer  term 
Tss .  The  reason  is  that  the  dominant  terms  in  Tss ,  e.g.,  the  incompressible  contributions  (Tfs  and  T^s), 
are  independent  of  Mt  (Figure  6).  On  the  other  hand,  the  “compressible”  contributions  (which  are  smaller) 
depend  strongly  on  the  values  of  the  turbulent  Mach  number.  Figure  7  illustrates  the  dependences  of  T§s  and 
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Fig.  11. 

Tfs  on  the  turbulent  Mach  number.  We  found  that  all  the  spectra  of  Tjf5  can  be  collapsed  by  dividing  M f 
(Figure  8).  Therefore,  we  conclude  that  T$s  is  proportional  to  M2,  a  result  that  can  be  found  analytically. 
Since  the  dominant  contributions  in  Tss  are  insensitive  to  the  variations  of  turbulent  Mach  number,  Tss  is 
not  affected  by  the  compressibility. 

3.2.  Compressible  Transfer  Term.  In  this  subsection,  we  will  study  the  transfer  term  Tcc,  and 
its  individual  contributors.  This  compressible  transfer  term  appears  in  the  transport  equation  of  the  com¬ 
pressible  auto-correlation.  As  shown  in  Figure  9,  the  compressible  transfer  is  positive  for  all  spectral  space. 
Hence,  it  is  a  term  that  is  responsible  for  the  production  of  compressible  energy. 

The  different  contributors  to  Tcc  are  plotted  in  Figure  10.  It  is  clear  that  the  two  terms,  T§c  and  Tjpc, 
are  much  larger  than  the  others.  Another  term,  Txcc,  is  a  distant  third  in  size.  To  illustrate  the  relative  size  of 
the  smaller  terms,  we  replotted  these  contributors  in  Figure  11  at  the  enlarged  scale.  Because  two  dominant 
terms  have  similar  magnitudes  but  opposite  signs,  a  strong  cancelation  between  them  is  expected.  Indeed 
for  all  cases  considered,  the  term  T§c  is  always  positive,  whereas  the  other  term  7]pc  is  always  negative. 
In  fact,  the  cancelations  are  so  complete,  the  summation  of  these  two  terms  is  now  actually  negligible  in 
comparison  with  Tpc  (Figure  12).  The  physical  explanation  for  this  ‘almost  perfect’  cancelation  is  that 
T%c  and  T«pc  are  the  terms  that  take  into  account  the  interactions  between  slowly  varying  incompressible 
modes  and  two  compressible  modes  (namely  the  interactions  between  acoustic  waves  and  a  solenoidal  field). 
This  interaction  results  in  the  production  of  acoustic  energy  on  the  same  wave-number  but  now  in  another 
propagation  direction.  For  an  isotropic  redistribution  of  acoustic  energy,  this  effect  of  reorientation  does  not 
affect  the  spectral  distribution  of  energy  and  leads  to  a  zero  net  balance.  Consequently,  the  most  important 
term  in  the  transport  equation  of  Ecc  is  T1cc.  Figure  13  demonstrates  that  this  term  is  also  much  larger 
than  all  other  contributors  (T3CC,  T%c  and  T^c). 
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Fig.  12. 

Comparing  TJ?S(K)  to  Tfc(K),  it  is  clear  that  they  have  similar  magnitude  but  with  an  opposite  sign 
(Figure  14).  These  two  terms  are  essentially  responsible  of  the  energy  exchanges  between  the  solenoidal  and 
compressible  parts.  Specifically,  Tfc  is  the  “input”  energy  term  on  the  compressible  mode  whereas  Tjf5  is 
the  “output”  term  in  the  equation  of  the  solenoidal  mode.  Based  on  these  results  obtained  in  spectral  space, 
we  conclude  that  there  is  a  local  transfer  of  energy  from  the  solenoidal  mode  to  the  compressible  mode.  This 
result  will  be  further  confirmed  by  our  analysis  in  the  second  part  of  the  paper. 

The  total  compressible  transfer  term  Tcc  is  dependent  on  the  compressibility  (see  Figure  15).  As 
expected,  its  magnitude  increases  with  the  increasing  of  turbulent  Mach  number.  We  note  that  there  is 
a  shift  of  the  peak  spectrum  towards  the  large  K .  Based  on  the  properties  of  T^S{K),  we  expect  that 
its  compressible  counterpart  Tfc(K)  should  also  scale  as  M2.  Since  the  spectra  of  divided  by  M2 
collapses  (Figure  16),  this  term  can  not  be  responsible  for  the  peak  shift  of  Tcc .  Figure  17  showed  that 
the  ‘almost  perfect’  cancelation  between  terms  T£c(K)  and  T^C{K)  for  all  turbulent  Mach  numbers  is 
considered.  Figure  18  illustrates  how  the  term  T^C(K)  depends  on  Mt .  Figure  19  shows  the  dependence 
of  T^C(K)  and  T§G(K)  on  the  values  of  turbulent  Mach  number.  Although  these  two  terms  maintain 
opposite  signs  for  all  spectral  space,  the  magnitudes  and  shape  of  T^C(K)  and  T§C(K)  are  clearly  different. 
The  ‘imperfect’  cancelation  between  these  two  terms  leads  to  a  cascade  type  of  compressible  energy  transfer 
(Figure  20).  These  two  terms  involve  ( Ecc )2  and  are  important  contributors  at  high  turbulent  Mach 
numbers.  The  interaction  among  the  compressible  mode  begins  to  have  influence.  This  cascade  mechanism 
will  be  investigated  in  the  next  section. 
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Fig.  13. 

4.  Triadic  Interactions.  The  most  fundamental  building  block  of  the  energy  transfer  process  is  the 
triadic  interactions.  Specifically,  we  are  interested  in  the  energy  transfer  for  a  given  mode  K  due  to  its 
interactions  with  all  the  pairs  of  modes  P  and  Q  =  K  -  P  that  form  a  triangle  with  K.  For  this  reason,  we 
introduce  the  triadic  energy  transfer  function,  T(K,  P,Q),  according  to 

(27)  TSS(K)=  J2  Tss (K,  P,Q),and 

P,Q= |K-P| 

(28)  TCC{K)  —  Tcc{K,P,Q). 

F,Q=|K-P| 

Here  T(K ,  P,  Q )  is  defined  as  energy  transfer  to  K  due  to  triads  with  one  leg  in  Q  and  the  other  in  P.  The 
average  procedure  is  performed  over  spherical  shell  since  the  turbulence  is  isotropic. 

An  examination  of  the  purely  incompressible  contributors  (TgS(K,  P,  Q))  reproduces  the  results  of  in¬ 
compressible  turbulence  (Domaradzki  and  Rogallo  (1990),  Yeung  et  al.  (1991,  1995),  Ohkitani  and  Kida 
(1992),  Zhou  (1993a-b)  and  Zhou  et  al.  (1996))  and  again  indicates  that  the  purely  solenoidal  triadic  energy 
transfer  is  not  affected  by  compressible  effects.  The  triadic  solenoidal  transfer  TSS(K ,  P,  Q)  (with  the  com¬ 
pressible  terms  included)  is  essentially  the  same  for  a  wide  range  of  turbulent  Mach  number  values  (Figure 
21).  Although  Figure  21  is  only  for  P  =  512  and  Q  —  128,  we  have  examined  other  values  of  P  and  Q  and 
found  that  our  conclusion  does  not  change  with  turbulent  Mach  number.  As  a  result,  the  compressibility 
has  very  little  influence  on  the  solenoidal  triadic  interactions.  Figure  22  is  a  typical  plot  showing  how  the 
structure  of  TSS(K,  P,Q)  changes  with  various  Q  values  (P  =  512,  Mt  =  10“2).  Again,  this  result  is  the 
same  as  that  of  incompressible  turbulence. 
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Figure  14:TSS5  andT^l 


Fig.  14. 


We  have  found  from  the  previous  section  that  particular  attention  should  be  paid  to  the  T§s  term  since 
it  is  the  term  that  is  responsible  for  ‘output’  energy  from  solenoidal  to  compressible  mode.  As  the  Mach 
number  increases,  the  magnitude  of  the  ‘output’  energy  increases  but  the  basic  structure  remains  (Figure 
23).  Since  this  term  represents  an  energy  output  at  a  localized  spectral  region,  we  refer  to  it  as  the  radiative 
(emission)  triadic  energy  transfer.  Figure  24  illustrates  this  type  of  interaction  from  various  Q  values  at 
the  given  value  of  P  (P  =  512)  for  Mt  =  10~2.  It  is  clear  that  the  triadic  interaction  of  this  term  is  quite 
different  from  those  of  purely  incompressible  terms. 

We  now  turn  our  attention  to  the  triadic  interactions  in  compressible  energy  transfer,  TCC(K\  P,Q). 
In  Figure  25,  we  present  TCC(K,  P,  Q)  for  various  Q  values  when  P  is  in  the  inertial  range  (P  =  512). 
For  this  low  Mach  number  10“2,  we  observe  that  the  structures  of  TCC{K,P,Q)  are  rather  similar  for 
differing  Q  values.  All  of  them  show  the  radiative  (absorption)  type  of  energy  transfer.  Figure  26  shows 
that  T^C(K,  P,  Q)  is  the  dominant  contributor  to  the  compressible  triadic  energy  transfer  function.  In 
fact,  comparing  Figures  24  and  26,  we  found  that  the  absorption  types  of  triadic  energy  transfer  functions 
have  the  same  magnitude  but  opposite  sign  as  those  of  the  emission  type  (T§s  term).  Figure  27  further 
demonstrates  that  the  triadic  interactions  T^C(K,  P,  Q)  &  —  T$S(K,  P,  Q)  for  all  turbulent  Mach  numbers 
are  under  consideration.  As  a  result,  we  conclude  that  all  compressible  energy  has  been  transferred  locally 
(in  spectral  space)  from  the  solenoidal  component. 

From  the  previous  section,  we  have  found  that  the  sum  of  T^C(K,  P,  Q)  +T^C(A",  P,  Q )  is  an  important 
contributor  to  the  compressible  energy  transfer.  At  low  Mach  number  ( Mt  =  10-2),  terms  T$C(K7  P,  Q)  + 
TqC{K,  P,  Q )  are  small  and  their  triadic  energy  transfer  terms  show  only  a  very  weak  energy  cascade  (Figure 
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Figure  15: Compressible  transfer  term  at  different  M( 


Fic?.  15. 


28).  However,  this  situation  changes  rapidly  as  the  Mach  number  increases.  Indeed,  the  compressible  energy 
cascade  can  be  seen  in  Figures  29  and  30  where  T^c  (K,  P,  Q )  +  T$C(K,  P,  Q )  are  plotted  for  several  higher 
Mach  numbers.  From  this  analysis,  we  conclude  that  at  high  Mach  number  the  cascade  of  compressible  tur¬ 
bulence  is  a  direct  result  of  the  fact  Tfc(K ,  P,  Q)  4-  Tfc(KJ  P,  Q)  >  Tfc(K,  P,  Q ).  To  further  demonstrate 
this  point,  we  plot  the  total  compressible  energy  transfer  term,  TCC(K ,  P,  Q),  at  different  Mach  numbers.  It 
is  clear  that  TCC(K ,  P,  Q )  changes  its  characteristic  features  from  radiative  to  cascade  as  the  turbulent  Mach 
number  increases  (Figures  31  and  32).  This  is  a  result  that  can  not  be  observed  from  the  total  compressible 
energy  transfer  function  TCC(K). 


14 


-  f3C+7f(Mt=10  2) 
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Figure  29:Triadic  T^+T^0  at  different  low  Mt 


Fig.  30. 
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Figure  30:Triadic  TC2C+1C6C  at  different  Mt 
Fig.  31. 
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-  fc{Mt=\Q~2) 

- fC(Mt= 10"1) 


101  102  103  104 
Figure  31:Triadic  T°c  at  different  low  Mt 


Fig.  32. 
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Figure  32:Triadic  1*~C  at  different  high  Mt 


Fig.  33. 

5.  Conclusion.  We  have  investigated  the  energy  transfer  process  of  compressible  turbulence  using  a 
two-point  closure  theory,  Eddy-Damped-Quasi-Normal-Markovian  (EDQNM)  closure,  as  well  as  the  Helmholtz 
decomposition,  a  method  that  separates  the  compressible  and  solenoidal  modes.  We  focused  on  the  following 
issues:  (1)  What  is  the  mechanism  of  energy  exchange  between  the  solenoidal  and  compressible  modes,  and 
(2)  Is  there  an  energy  cascade  in  the  compressible  energy  transfer  process?  We  found  that  the  compressible 
energy  is  transferred  locally  from  the  solenoidal  part  to  the  compressible  part.  We  also  found  that  there  is 
an  energy  cascade  of  the  compressible  mode  for  high  turbulent  Mach  number  ( Mt  >  0.5).  Since  we  assume 
that  the  compressibility  is  weak,  the  magnitude  of  the  compressible  (radiative  or  cascade)  transfer  is  much 
smaller  than  that  of  solenoidal  cascade.  We  confirmed  these  results  by  studying  the  triadic  energy  transfer 
function,  a  most  fundamental  aspect  of  the  energy  transfer. 
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